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Abstract. We study the possible spin-1 charmonium-like states by using QCD sum rule 
approach. We calculate the two-point correlation functions for all the local form tetraquark 
interpolating currents with J FL = 1 , L + , 1 ++ and l + ~ and extract the masses of the 
tetraquark charmonium-like states. The mass of the 1 qcqc state is 4.6 ~ 4.7 GeV, which 
implies a possible tetraquark interpretation for 7(4660) meson. The masses for both the 
1 ++ qcqc and scsc states are 4.0 ~ 4.2 GeV, which is slightly above the mass of X(3872). 
For the 1 -+ and l + ~ qcqc states, the extracted masses are 4.5 ~ 4.7 GeV and 4.0 ~ 4.2 
GeV respectively. 
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Qu 1 Introduction 



The underlying structures of the so-called X, Y, Z states are not understood well I1I2I31 . Some of them 
do not fit in the conventional quark model easily and are considered as the candidates of the exotic 
states beyond the quark model, such as the molecular states, tetraquark states, the charmonium hy- 
brid mesons, baryonium states and so on. X(3872) is the best studied charmonium-like state since its 
discovery by the Belle Collaboration (4). Although the analysis of angular distributions favors the as- 
signment/ PC = 1 ++ El?], the 2~ + possibility is not ruled out [6|. The mass and decay mode of X(3872) 
are very different from that of the 2^P\ cc state. Up to now, the possible interpretations of X(3872) 
include the molecular state I7I8I9I10I . tetraquark state 11 11121 . cusp Ifl3l and hybrid charmonium 1141 . 
7(4260), F(4360) and 7(4660) are the Y(J PC = 1") family states discovered in the initial state radia- 
tion (ISR) process. These new states lie above the open charm threshold. However, the Y — » 
decay modes have not been observed yet 03], which are predicted to be the dominant decay modes of 
the charmonium above the open charm threshold in the potential model. In Refs. 1161171 . the authors 
studied the 1 charmonium-like Y mesons using the QCD sum rule approach. Maiani et al. tried to 
assign 7(4260) as the scsc tetraquark in a P-wave state 11181 . 7(4260) was also interpreted as the in- 
teresting charmonium hybrid state II 191201 . 7(4660) was considered as a (/r(25)/o(980) bound state in 
Ref. (211. Recently there have been some efforts on the 1~ + charmonium-like exotic states. For exam- 
ple, the structure of X(4350) was studied using a D*D* Q current with J pc — 1~ + [22 1 . Moreover, the 
newly observed state 7(4140) was argued as a 1~ + exotic charmonium hybrid state |23l . 

In this contribution, we would like to report the systematic study of the tetraquark charmonium-like 
states with J PC = 1 , 1~ + , 1 ++ and l + ~. After constructing all the tetraquark interpolating currents 
with definite quantum numbers, we investigate the two-point correlation functions and extract the 
masses of the charmonium-like states with QCD sum rule. We study both the qQqQ and sQsQ systems 
where Q — c,b. We also discuss the possible decay modes and experimental search of the charmonium- 
like states. 
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2 INDEPENDENT CURRENTS 



Firstly, we construct all the local form diquark-antidiquark type of interpolating currents in a sys- 
tematic way. By considering the Lorentz structures, the charge conjugation properties and the color 
symmetries of the tetraquark operators, we arrive the following tetraquark interpolating currents with 
J p c = 1 +, 1 — , 1++ and 

- The interpolating currents with J pc — 1~ + and 1 are: 

J in = qlCysCbiqay^iCcl + q b y^y 5 Ccl) + q T a Cy ll y^c b {q a y' i Cc T b + q b ysCc T a ) , 

hn = q T a Cy v c b {q a a- ilv Cc T b - qbCr^Ccl) ± q T a Ccr tlv c b {q a y y Cc T b - q b y v Cc T a ) , 

hp = qlCysCbiqaJi/YsCcl - qbJ/YsCcl) + qlCy M y 5 c b (q a y 5 Ccl - q b y 5 Cc T a ) , 

Jfy = qlCy v c b (q a o- flv Ccl + q b o-^ v Ccl) ± q T a Ccr llv c b (q a y v Cc 1 b + q b y v Cc T a ) , 

hp = llCc b (q a y M Ccl + q b y M Ccl) + qlCy^CbiqaPcl + q b Cc T a ) , (1) 

J 6fl = qlCy v y 5 c b (q a o- flv y 5 Ccl + qhCr^ysCcl) + q T a C(T llv y i c b {q a y v y i Cc T b + q b y v ysCc T a ) , 

hp = qlCc b (q a y M Ccl - q^Cc^) + q T a Cyy.c b {q a Cc T b - q b Cc T a ) , 

hp = qlCy v y 5 c b (q a o- flv y 5 Ccl - qhO-^y^Ccl) + q T a C(T ilv y i c b {q a y v ysCc l: b - q b y v y 5 Ccl) . 

where "+" corresponds to J PC = "-" corresponds to J pc = 1 . 

- The interpolating currents with J pc — 1 ++ and l + ~ are: 

J in = q T a Cc b (q a ypy$Cc T b + qbJpysCcl) + qlCy^ysChiqaCcl + q b Cc T a ) , 

hp = qlCc b (q a y M y 5 Cc T b - q b y M y 5 Ccl) + q T a Cypy$c b {q a Cc T b - q b Cc T a ) , 

hp = qlCy 5 c b (q a y M Ccl + q b y M Ccl) + q T a CypC b {q a ysCc T b + q b y s Cc T a ) , 

hp = qlCy 5 c b (q a y M Ccl - q b y M Ccl) + q T a CypC b (q a y 5 Cc T b - q b y s Cc T a ) , (2) 

J 5fI = q T a Cy v c b {q a crn V yiCc T b + q b o-p V y 5 Ccl) ± q T a Ca-n V yiC b {q a y v Cc T b + q b y v Cc T a ) , 

hp = qlCy v c b (q a o- flv y 5 Cc T b - q b CTp V y 5 Ccl) + q T a Ccr ll vy : iC b {q a y v Cc 1 b - q b y y Cc T a ) , 

J lfl = q T a Cy v ysc b {q a o-n V Cc T b + q b o-p V Ccl) ± q T a C(XnvC b {q a y v ysCc 1 b + q b y v y 5 Ccl) , 

hp = qJ,Cy v y 5 c b (q a o- flv Cc T b - qhCr^Ccl) + qlCa llv c b (q a y v y5Ccl - q b y y ysCc T a ) . 

where "+" corresponds to J PC - 1 ++ , "-" corresponds to J PC - 

The subscripts a and b are the color indices, q denotes u or d quark. It is understood that all the currents 
in Eqs. (Hi-© should contain (ucuc + dcdc) in order to have definite isospin and G-parity. The details 
about the current construction could be found in Ref . Il24l . 



In the past several decades, QCD sum rule has been widely used to study the hadron structures and 
proven to be a very powerful non-perturbative method 0251261 . We consider the two-point correlation 
function: 



3 SPECTRAL DENSITY 




(3) 



where J M is a interpolating current for the tetraquark states. II\{q 2 ) is related to the vector meson while 
IJo(q 2 ) is the scalar current polarization function. The correlation function n^viq 2 ) can be calculated 
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in the operator product expansion (OPE) using perturbative QCD augmented with non-perturbative 
quark and gluon condensates to describe the large distance physics. At the hadron level, the correlation 
function is expressed by the dispersion relation with a spectral function: 



JAml s-q- -le 



In approximation of the infinitely narrow widths of resonances, the spectral function can be expressed 
as: 



p(s) = Y j d(s-ml)(0\r ] \n)(n\r, + \0) 



= f 2 6(s - m\) + (5) 

where represents the higher states contribution. 

The theoretical basis of the QCD sum rule approach is the assumption of the quark-hadron duality, 
which ensures the equivalence of the correlation functions obtained at the hadron level and the quark- 
gluon level. After performing the Borel transformation to the correlation functions, we can extract the 
mass of the state X: 

m x = Tin 1 77772 ~~- (6) 



fif 1 — L 



where so is the continuum threshold and Mb is the Borel parameter. We performed the QCD sum rule 
analysis for all the tetraquark currents in Eqs. JO-©. The results of OPE can be found in Ref. 11241 . 



4 QCD Sum Rule Analysis 

We use the following parameter values of the quark masses and various condensates 1271281291 for 
the numerical analysis: m 9 (2GeV) = (4.0 + 0.7) MeV,m J (2GeV) = (lOl^)MeV, m c (m c ) = (1.23 ± 
0.09)GeV, m b (m b ) = (4.20 + 0.07)GeV, (qq) = -(0.23 + 0.03) 3 GeV 3 , (qg s crGq) = -Ml(qq),M 2 = 
(0.8 + 0.2)GeV 2 , (ss)/(qq) = 0.8 + 0.2, (g 2 ,GG) = 0.88GeV 4 . There are two important parameters 
in QCD sum rule analysis: the threshold parameter so and the Borel mass Mb- The stability of QCD 
sum rule requires a suitable working region of sq and Mb- Since the exponential weight function 
in Eq. |6l the higher state contribution is naturally suppressed for small value of Mb- However, the 
OPE convergence would become worse if Mb was too small. In our analysis, we choose the value of 
so around which the variation of the extracted mass mx with M 2 B is minimum. The working region 
of the Borel mass is determined by the convergence of the OPE series and the pole contribution. The 
requirement of the convergence of the OPE series leads to the lower bound M 2 . of the Borel parameter 
while the constraint of the pole contribution yields the upper bound of M 2 B . 



4.1 Vector charmonium-like systems 

For the interpolating currents with J PC — 1~ + and 1 , we keep the m q and m s related terms in the 
spectral densities. These terms give important corrections to the OPE series and are useful to enhance 
the stability of the sum rule. For the qcqc systems, the absolute value of the four quark condensate 
{qq) 2 is bigger than other condensates in the region of M 2 < 3.1 GeV 2 . It is the dominant power 
contribution to the correlation function in this region. Especially for the currents with J PC = 1 , the 
quark condensate (qq) is proportional to the light quark mass m q and vanishes if we take m q = 0. 
However, it is proportional to the strange quark mass m s and larger than (ss) 2 for the scsc system 
since m s » m q . This is the main difference between the qcqc and scsc systems. The similar situation 
exits in 1~ + charmonium-like systems. 
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After careful study of the OPE convergence and the pole contribution, we find the suitable working 
region of the Borel parameter for each vector charmonium-like current. The threshold value of so is 
also fixed around which the variation of the extracted mass nix with M 2 B is minimum. In Fig. Q] and 
Fig-El we show the variation of nix with the threshold value so and Borel parameter M 2 B for the current 
Ji^ with J pc - 1 in qcqc and scsc systems, respectively. One notes that they are very similar with 
each other except the chosen so mentioned above. The extracted mass of the qcqc state is 4.64 GeV, 
which is consistent with the mass of the meson T(4660). One may wonder whether T(4660) could be 
a tetraquark state. The extracted mass of the scsc state is 4.92 GeV, which is about 0.28 GeV higher 
than that of the qcqc state. 
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Fig. 1. The variation of m x with ,v (a) and M~ B (b) corresponding to the current J lfl for the 1 qcqc system. 





5.5 




5.0 


> 


4.5 


O 








1 


4.0 




3.5 




3.0 







%> 
t 

K 


















Af|=3.0 

Mj5=3.2 

Mb =3 .4 




Hi'. 

i\Y. 




M*=3.6 



5.5 
5.0 



% 4.5 

p 



I 4 -° 
3.5 

3.0 



ll \ . 

\\ \ \ 
1 * 1 \ 
' \ \ '• . 
\\ \ 






Mb =2. 8 




Mg=3.3 




«a=3.8 








Af |=4.3 











5.2 I 

5.0 

4.8 - 

a 

„ 4.6 ,, 0= 5.2 2 

S 4.4- I0=5 .4 ! 

4.2 - ,, 0= 5.6 2 

4.0 - 



12 15 18 21 24 27 30 33 2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 

s [GeV 2 ] Ml [GeV 2 ] 

(a) (b) 
Fig. 2. The variation of m x with ^o( a ) an d M 2 B (Jo) corresponding to the current for the 1 scsc system. 

Performing the QCD sum rule analysis, we show the Borel window, the threshold value, the ex- 
tracted mass and the pole contribution corresponding to the tetraquark currents with J PC = 1 in 
Table [T] The results of the 1~ + system are listed in Table [2] We only present the numerical results for 
the currents which lead to the stable mass sum rules in the working range of the Borel parameter. For 
example, only the currents Ji^, and J-/^ with J pc = 1 in qcqc systems have the reliable mass sum 
rules. For J^^, J^, Js^, and J^, the stability is so bad that the extracted mass nix grows monoton- 
ically with the threshold value so and the Borel parameter Mb- These currents may couple to the 1 
states very weakly, leading to the above unstable mass sum rules. We also study the bottomonium-like 
analogues by replacing m c with nib in the correlation functions and repeating the same analysis pro- 
cedures done above. The numerical results of the qbqb and sbsb systems are collected in Table Q] and 
Table|2] 



Hadron Nuclear Physics (HNP) 201 1 





Currents 


so(GeV-) 


[M 2 . ,M 2 ](GeV 2 ) 
mm max JV * 


ra x (GeV) 


PC(%) 




j, 

J i/i 


j.yj 




A f\A + OQ 
t.Ut- it \j.\jy 




qcqc system 


J An 


5.0 2 


2.9 ~ 3.6 


4.61 +0.10 


46.4 




hp 


5.2 2 


2.9-4.1 


4.74 + 0.10 


47.3 




hp 


5.4 2 


2.8-4.5 


4.92 + 0.10 


50.3 




hp 


5.0 2 


2.8 ~ 3.5 


4.64 ± 0.09 


48.6 


scsc system 


hp 


4.9 2 


2.8 ~ 3.4 


4.52 + 0.10 


45.6 




hp 


5.4 2 


2.8-4.5 


4.88+0.10 


51.7 




hp 


5.3 2 


2.8-4.3 


4.86 + 0.10 


46.0 




hp 


4.8 2 


2.8 - 3.1 


4.48+0.10 


43.2 


qbqb system 


hp 


11.0 2 


7.2 - 8.5 


10.51+0.10 


45.8 




hp 


11.0 2 


7.2 - 8.3 


10.60 + 0.10 


47.0 




hji 


11.0 2 


7.2 - 8.4 


10.55 + 0.11 


43.6 


sbsb system 


hp 


11.0 2 


7.2 - 8.4 


10.55 + 0.10 


43.7 




J Ail 


11.0 2 


7.2 - 8.4 


10.53 + 0.11 


44.3 




hji 


11.0 2 


7.2 - 8.2 


10.62 + 0.10 


42.0 




hp 


11.0 2 


7.2 - 8.4 


10.53 + 0.10 


44.1 



Table 1. The threshold value, Borel window, mass and pole contribution corresponding to the currents with 
J FC = 1 in the qcqc, scsc, qbqb and sbsb systems. 





Currents 


so(GeV 2 ) 


[^in^fnaxKGeV 2 ) 


m x (GeV) 


PC(%) 




hp 


5.1 2 


2.9 - 


3.9 




4.67 ±0.10 


50.2 


qcqc system 


hp 


5.2 2 


2.9 - 


4.2 




4.77+0.10 


47.4 




Jgfi 


4.9 2 


2.9 - 


3.4 




4.53+0.10 


46.3 




Jip 


5.0 2 


2.9 - 


3.4 




4.67+0.10 


44.3 




hp 


5.0 2 


2.9 - 


3.4 




4.65 ± 0.09 


45.6 




hp 


4.9 2 


2.9 - 


3.3 




4.54 + 0.10 


44.4 






5.1 2 


2.9 - 


3.7 




4.72 ± 0.09 


44.8 


scsc system 


hp 


5.0 2 


2.9 - 


3.6 




4.62 + 0.10 


42.8 




hp 


5.3 2 


2.9 - 


4.3 




4.84 + 0.10 


47.3 




hp 


5.3 2 


2.9 - 


4.3 




4.87 + 0.10 


46.2 




hp 


5.2 2 


2.9 - 


4.1 




4.77 + 0.10 


44.1 




hp 


n.o 2 


7.2 - 


8.6 




10.53 + 0.11 


44.2 


qbqb system 


hp 


ll.O 2 


7.2 - 


8.6 




10.53 + 0.10 


44.1 




hp 


ll.O 2 


7.2 - 


8.6 




10.49 + 0.11 


44.7 




hp 


ll.O 2 


7.2 - 


8.1 




10.62 + 0.10 


41.2 




hp 


ll.O 2 


7.2 - 


8.4 




10.56 + 0.10 


43.8 


qbqb system 


hp 


ll.O 2 


7.2 - 


8.3 




10.63 + 0.10 


42.4 




hp 


ll.O 2 


7.2 - 


8.3 




10.62 ± 0.09 


42.5 




hp 


ll.O 2 


7.2 - 


8.3 




10.59 + 0.10 


43.1 



Table 2. The threshold value, Borel window, mass and pole contribution corresponding to the currents with 
J PC = 1~ + in the qcqc, scsc, qbqb and sbsb systems. 

4.2 Axial-vector charmonium-like systems 

In this channel, the QCD sum rule analysis shows that the quark condensate (qq) is the dominant 
correction to the correlation function for all the currents. Although the OPE convergence becomes 
worse than that in the vector channel, the currents J^, J^, J-j^, have better OPE convergence than 
that of Ji^Ji^J^Jap- For the interpolating currents and with J PC = 1 ++ , we obtain the 
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working region of the Borel parameter in 3.0 < M\ < 3.4 GeV 2 for both the qcqc and scsc systems. 
The extracted mass is about mx = 4.0 ~ 4.2 GeV, which is slightly above the mass of X(3872). The 
qbqb and sbsb systems can be studied conveniently by replacement of the parameters, including the 
quark masses and the various condensates. The numerical results are listed in Table [3] for the 1 ++ 
systems and Table|4]for the l + ~ systems. 





Currents 


so(GeV-) 


[^ min ^max](GeV 2 ) 


m x (GeV) 


PC(%) 


qcqc system 


Jin 


4.6 2 


3.0 ~ 


3.4 


4.19 + 0.10 


47.3 




J Ail 


4.5 2 


3.0 ~ 


3.3 


4.03+0.11 


46.8 


scsc system 




4.6 2 


3.0 ~ 


3.4 


4.22 + 0.10 


45.7 




J Ail 


4.5 2 


3.0 ~ 


3.3 


4.07 + 0.10 


44.4 




J in 


10.9 2 


8.5 ~ 


9.5 


10.32 ± 0.09 


47.0 


qbqb system 


J Ail 


10.8 2 


8.5 ~ 


9.2 


10.22 + 0.11 


44.6 




■hn 


10.7 2 


7.8 ~ 


8.4 


10.14 + 0.10 


44.8 




Hn 


10.7 2 


7.8 ~ 


8.4 


10.14 + 0.09 


44.8 




hii 


10.9 2 


8.5 ~ 


9.5 


10.34 ± 0.09 


46.1 


sbsb system 


J Ail 


10.8 2 


8.5 ~ 


9.1 


10.25 + 0.10 


43.3 




Jiii 


10.8 2 


7.5 ~ 


8.6 


10.24 + 0.11 


47.1 




J in 


10.8 2 


7.5 ~ 


8.6 


10.24 + 0.10 


47.1 



Table 3. The threshold value, Borel window, mass and pole contribution corresponding to the currents with 
J FC = 1 ++ in the qcqc, scsc, qbqb and sbsb systems. 





Currents 


so(GeV 2 ) 


[^ min ^fna X ](GeV 2 ) 


m x (GeV) 


PC(%) 




Jin 


4.6 2 


3.0 ~ 3.4 




4.16 + 0.10 


46.2 


qcqc system 


J Ail 


4.5 2 


3.0 ~ 3.3 




4.02 ± 0.09 


44.6 




■hi, 


4.5 2 


3.0 ~ 3.4 




4.00 + 0.11 


46.0 




Jbn 


4.6 2 


3.0 ~ 3.4 




4.14 + 0.09 


47.0 




Jiii 


4.7 2 


3.0 ~ 3.6 




4.24 + 0.10 


49.6 


scsc system 


J Ail 


4.6 2 


3.0 ~ 3.5 




4.12 + 0.11 


47.3 




Jsn 


4.5 2 


3.0 ~ 3.3 




4.03+0.11 


44.2 




J6li 


4.6 2 


3.0 ~ 3.4 




4.16 + 0.11 


46.0 




Jin 


10.6 2 


7.5 ~ 8.5 




10.08 + 0.10 


45.9 


qbqb system 


J Ail 


10.6 2 


7.5 ~ 8.5 




10.07 + 0.10 


46.2 




Jsn 


10.6 2 


7.5 ~ 8.4 




10.05 + 0.10 


45.3 




J6pi 


10.7 2 


7.5 ~ 8.7 




10.15 + 0.10 


47.6 




J in 


10.6 2 


7.5 ~ 8.3 




10.11+0.10 


43.8 


sbsb system 


J Ail 


10.6 2 


7.5 ~ 8.4 




10.10 + 0.10 


44.1 




J 5/1 


10.6 2 


7.5 ~ 8.3 




10.08 + 0.10 


43.7 




J 6/i 


10.7 2 


7.5 ~ 8.5 




10.18 + 0.10 


46.5 



Table 4. The threshold value, Borel window, mass and pole contribution corresponding to the currents with 
J pc = l + ~ in the qcqc, scsc, qbqb and sbsb systems. 



5 Conclusion 

We have performed the QCD sum rule analysis with tetraquark charmonium-like currents in vector 
and axial-vector channels. The two-point correlation functions and the spectral densities for all the 
interpolating currents have been evaluated at the quark-hadron level. The numerical analysis shows 
that the four quark condensate (qq) 2 is the dominant power contribution to the OPE series for all 
the vector channel currents. In the situation of the axial-vector channel currents, however, the most 
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important corrections are the quark condensate (qq). These properties of the spectral densities lead to 
a better OPE convergence for the currents in the vector channel than that in the axial-vector channel. 
The m s related terms in the OPE series of scsc systems lead to more stable mass sum rules than that 
of the qcqc systems. In the working range of the Borel parameter, only the currents Ji M , and Jj M 
with J PC = 1 display stable QCD sum rules in the qcqc system. The extracted mass is around 
4.6 ~ 4.7 GeV from these currents, which is consistent with the mass of the meson F(4660). This 
result implies a possible tetraquark interpretation for T(4660). In the scsc system, all currents except 
J511, J(,n have stable QCD sum rules and the extracted mass is about 4.6 ~ 4.9 GeV. The Borel window 
for the currents in the axial-vector channel is very small because of the bad OPE convergence. For the 
currents with J pc = 1 ++ in the qcqc system, only and have reliable QCD sum rules. The same 
situation occurs in the scsc system. The extracted masses are about 4.0 ~ 4.2 GeV, which is 0.1 ~ 0.3 
GeV higher than the mass of X(3872). 

The possible decay modes of these charmonium-like states are also studied by considering the 
conservation of the angular momentum, P-parity, C-parity, isospin and G-parity ||24|| . 
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